Lesion Border Detection in Dermoscopy Images 
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Abstract 

Background: Dermoscopy is one of the major imaging modalities used in the diagnosis of melanoma and other pigmented 
skin lesions. Due to the difficulty and subjectivity of human interpretation, computerized analysis of dermoscopy images has 
become an important research area. One of the most important steps in dermoscopy image analysis is the automated detection 
of lesion borders. Methods: In this article, we present a systematic overview of the recent border detection methods in 
the literature paying particular attention to computational issues and evaluation aspects. Conclusion: Common problems 
with the existing approaches include the acquisition, size, and diagnostic distribution of the test image set, the evaluation 
of the results, and the inadequate description of the employed methods. Border determination by dermatologists appears to 
depend upon higher-level knowledge, therefore it is likely that the incorporation of domain knowledge in automated methods 
will enable them to perform better, especially in sets of images with a variety of diagnoses. 



1 Introduction 

Invasive and in-situ malignant melanoma together comprise one of the most rapidly increasing cancers in the world. Invasive 
melanoma alone has an estimated incidence of 62,480 and an estimated total of 8,420 deaths in the United States in 2008 [TJ. 
Early diagnosis is particularly important since melanoma can be cured with a simple excision if detected early. 

Dermoscopy, also known as cpilumincscence microscopy, is a non-invasive skin imaging technique that uses optical magnifi- 
cation and either liquid immersion and low angle-of-incidence lighting or cross-polarized lighting, making subsurface structures 
more easily visible when compared to conventional clinical images [2] . Dermoscopy allows the identification of dozens of mor- 
phological features such as pigment networks, dots/globules, streaks, blue- white areas, and blotches [3]. This reduces screening 
errors, and provides greater differentiation between difficult lesions such as pigmented Spitz nevi and small, clinically equivocal 
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lesions [4]. However, it has been demonstrated that dermoscopy may actually lower the diagnostic accuracy in the hands 
of inexperienced dermatologists [5]. Therefore, in order to minimize the diagnostic errors that result from the difficulty and 
subjectivity of visual interpretation, the development of computerized image analysis techniques is of paramount importance [3] . 

Automated border detection is often the first step in the automated or semi-automated analysis of dermoscopy images [7] . 
It is crucial for the image analysis for two main reasons. First, the border structure provides important information for accurate 
diagnosis, as many clinical features such as asymmetry, border irregularity, and abrupt border cutoff are calculated directly 
from the border. Second, the extraction of other important clinical features such as an atypical pigment network, globules, 
and blue-white areas, critically depends on the accuracy of border detection. Automated border detection is a challenging 



task due to several reasons: (i) low contrast between the lesion and the surrounding skin (Fig. 1(a)), (ii) irregular (Fig. 1(b)) 



and fuzzy lesion borders (Fig. 1(c)), (hi) artifacts and intrinsic cutaneous features such as black frames, skin lines, blood 



vessels (Fig. 1(d)), hairs (Fig. 1(e)), and air bubbles (Fig. 1(f)), (iv) variegated coloring inside the lesion (Fig. 1(g)), and 



(v) fragmentation due to various reasons such as scar-like depigmentation (Fig. 1(h) ) 




(a) Low contrast 



(b) Irregular border 



(c) Fuzzy border 



(d) Blood vessels 




(f) Bubbles (g) Variegated coloring 

Figure 1: Problems with border detection 



(h) Fragmentation 



In the last decade, numerous methods have been developed for automated border detection in dermoscopy images. In 
this article, we present an overview of recent border detection methods and describe the preprocessing, segmentation, and 
postprocessing steps involved in each method. We also review performance evaluation issues and propose guidelines for studies 
in automated border detection. 



2 Preprocessing 

In this section, we describe the preprocessing steps that facilitate the border detection procedure namely, color space transfor- 
mation, contrast enhancement, and artifact removal. 

2.1 Color Space Transformation 

Dermoscopy images are commonly acquired using a digital camera with a dermoscope attachment. Due to the computational 
simplicity and convenience of scalar (single channel) processing, the resulting RGB (red-green-blue) color image is often 
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converted to a scalar image using one of the following methods: 

• Retaining only the blue channel (lesions are often more prominent in this channel). 

• Applying the luminance transformation, i.e. Luminance = 0.299 x Red + 0.587 x Green + 0.114 x Blue. 

• Applying the Karhunen-Loeve (KL) transformation [8] and retaining the channel with the highest variance. 

In applications where vector (multichannel) processing is desired, the RGB image can be used directly or it might be transformed 
into a different color space for various reasons including: (i) reducing the number of channels, (ii) decoupling luminance and 
chromaticity information, (iii) ensuring approximate perceptual uniformity, and (iv) achieving invariancc to different imaging 
conditions such as viewing direction, illumination intensity, and highlights. Common target color spaces in this case include 
CIEL*a*b*, CIEL*u*v*, KL, and HSI (Huc-Saturation-Intcnsity) [8]. Note that the use of the CIEL*a*b* and CIEL*u*v* 
color spaces requires careful calibration of the acquisition device, a step that seems to be frequently neglected in the literature. 

2.2 Contrast Enhancement 

As mentioned in Section [TJ one of the factors that complicate the detection of borders in dermoscopy images is insufficient 
contrast. Recently, Dclgado et al. |jJJ proposed a contrast enhancement method based on independent histogram pursuit (IHP). 
This algorithm linearly transforms the original RGB image to a decorrelated color space in which the lesion and the background 
skin are maximally separated. Border detection is then performed on these contrast-enhanced images using a simple clustering 
algorithm. 

2.3 Artifact Removal 

Dermoscopy images often contain artifacts such as such as black frames, ink markings, rulers, air bubbles, as well as intrinsic 
cutaneous features that can affect border detection such as blood vessels, hairs, and skin lines. These artifacts and extraneous 
elements complicate the border detection procedure, which results in loss of accuracy as well as an increase in computational 
time. The most straightforward way to remove these artifacts is to smooth the image using a general purpose filter such as 
the Gaussian (GF), median (MF), or anisotropic diffusion filters (ADF). Several issues should be considered while using these 
filters: 

• Scalar vs. vector processing: These filters are originally formulated for scalar images. For vector images one can apply 
a scalar filter to each channel independently and then combine the results, a strategy referred to as marginal filtering. 
Although fast, this scheme introduces color artifacts in the output. An alternative solution is to use filters that treat the 
pixels as vectors [TU] , 

• Mask size: The amount of smoothing is proportional to the mask size. However, excessively large masks result in the 
blurring of edges, which might reduce the border detection accuracy. Setting the mask size proportional to the image 
size seems to be a reasonable strategy [TTJ [T^] . 

• Computational time: For the GF and MF, algorithms that perform in constant time independent of the mask size have 
been developed [131 HI]. As for the ADF, the computational time depends on the mask size and the number of iterations. 
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An alternative strategy for artifact removal is to use specialized methods for each artifact type. For the removal of black frames, 
Celebi et al. [12] proposed an iterative algorithm based on the lightness component of the HSL (Huc-Saturation-Lightness) 
color space. In most cases, image smoothing effectively removes the skin lines and blood vessels. Hair removal received the 
most attention in the literature. Lee et al. [15) and Schmid |llj used mathematical morphology. Fleming et al. [BJ applied 
curvilinear structure detection with various constraints followed by gap filling. Recently, Zhou et al. [16] and Wighton et 
al. |17| proposed more sophisticated approaches based on inpainting. A method that can remove bubbles with bright edges 
was introduced in [BJ, where the authors utilized a morphological top-hat operator followed by a radial search procedure. 

3 Segmentation 

Segmentation refers to the partitioning of an image into disjoint regions that are homogeneous with respect to a chosen property 
such as luminance, color, texture, etc. |18| . Segmentation methods can be roughly classified into the following categories: 

• Histogram thresholding: These methods involve the determination of one or more histogram threshold values that separate 
the objects from the background. 

• Clustering: These methods involve the partitioning of a color (feature) space into homogeneous regions using unsupervised 
clustering algorithms. 

• Edge-based: These methods involve the detection of edges between the regions using edge operators. 

• Region-based: These methods involve the grouping of pixels into homogeneous regions using region merging, region 
splitting, or both. 

• Morphological: These methods involve the detection of object contours from predetermined seeds using the watershed 
transform. 

• Model-based: These methods involve the modeling of images as random fields whose parameters are determined using 
various optimization procedures. 

• Active contours (snakes and their variants): These methods involve the detection of object contours using curve evolution 
techniques. 

• Soft computing: These methods involve the classification of pixels using soft-computing techniques including neural 
networks, fuzzy logic, and evolutionary computation. 

Several issues should be considered when choosing a segmentation method: 

• Scalar vs. vector processing: Most segmentation methods are designed for scalar images. Although numerous vector 
image segmentation methods have been developed in the past decade, their usage is hindered by various factors including 
excessive computational time requirements and the difficulty of choosing an appropriate color space. 

• Automatic vs. semi-automatic: Some segmentation methods are completely automated, whereas others require human 
interaction. For example, active contour methods often require the manual delineation of the initial contour, whereas 
seeded region growing methods require the specification of the initial region seeds. 
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• Number of parameters: Most segmentation methods have several parameters whose values need to be determined a priori. 
In general, the more the number of parameters, the harder the model selection (determination of the optimal parameters). 

4 Postprocessing 

The result of the segmentation procedure is either a label image or a binary edge map. In order to obtain the lesion border, the 
segmentation output should be postprocesscd. The precise sequence of postprocessing operations depends on the particular 
choice of the segmentation method. However, certain operations seem to be useful in general. These include: 

• Region merging: Ideally, the segmentation procedure is expected to produce two regions: the lesion and the background 
skin. However, since these regions are rarely homogeneous, segmentation methods often partition them into multiple 
subregions. In order to obtain a single lesion object, subregions that are part of the lesion should first be identified and 
then merged. This can accomplished in several ways: 

— If the black frame of the image has already been removed, the background skin color can be estimated from the 
corners of the image and the subregions with similar color to the background skin can be eliminated, leaving only 
those subregions that are part of the lesion JT2l [19l [20] . 

— Various color and texture features can be extracted from each region and a classifier can be trained to determine 
which features effectively discriminate between the regions that are part of the lesion and those that are part of the 
background skin [5J. 

— The partitioning that maximizes the normalized texture gradient at the border and the total pairwise similarity 
among the regions within the lesion and those within the background skin can be determined by minimizing a cost 
function [16] , 

• Island removal: Islands (small isolated regions) in the label image can be eliminated using a binary area opening filter. 

• Border smoothing: Most segmentation methods produce regions with ragged borders. More natural borders can be 
obtained by a variety of operations including majority filtering |12j . morphological filtering [11[ [9]. and curve fitting. 

• Border expansion: In several studies, it was observed that the computer-detected borders were mostly contained within 
the dermatologist-determined borders. This is because the automated segmentation methods tend to find the sharpest 
pigment change, whereas the dermatologists choose the outmost detectable pigment. The discrepancy between the two 
borders can be reduced by expanding the computer-detected border using morphological filtering |12] , Euclidean distance 
transform |12j . or iterative region growing |21j . 

5 Evaluation 



Evaluation of the results seems to be one of the least explored aspects of the border detection task. As in the case of the more 
general image segmentation problem, there are two major evaluation methods: subjective and objective. The former involves 
the visual assessment of the border detection results by one or more dermatologists. Since there is no objective measure 
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of quality involved, this technique does not permit parameter tuning or comparisons among automated border detection 
methods. On the other hand, objective evaluation involves the quantification of the border detection errors using dermatologist- 
determined borders. In the rest of this discussion, we refer to the computer-detected borders as automatic borders and 
those determined by dermatologists as manual borders. Most of the quantitative error metrics are based on the concepts of 
true/false positive/negative given in Table [1] (here actual and detected pixels refer to a pixel in the ground-truth image and 
the corresponding pixel in the border detection output, respectively). These include: 

• XOR measure |32] = x 100% 

• Sensitivity = rpp~^pj\ j X 100% and Specificity = pp + 'pj\ | x 100% 

• Precision (Positive Predictive Value) = ^pp^pp x 100% and Recall = Sensitivity 

• Error probability = r pp + p " ]^p " p + r pj\ | x 100% 



Table 1: Definitions of true/false positive/negative 





Detected Pixel 


Lesion 


Background 


Actual 
Pixel 


Lesion 


True Pos. (TP) 


False Ncg. (FN) 


Background 


False Pos. (FP) 


True Neg. (TN) 



In a comprehensive study, Guillod et al. |23j demonstrated that a single dermatologist, even one who is experienced in 
dermoscopy, cannot be used as an absolute reference for evaluating border detection accuracy. In addition, they emphasized 
that manual borders are not precise, with inter-dermatologist borders and even borders determined by the same dermatologist 
at different times showing significant disagreement, so that a probabilistic model of the border is preferred to an absolute 
gold-standard model. Accordingly, they used fifteen sets of borders drawn by five dermatologists over a minimum period of one 
month. A probability image for each lesion was constructed by associating a misclassification probability p(i,j) = 1 — "^ J ' > 
with each pixel (N: number of observations, n(i, j): number of times pixel was selected as part of the lesion). For each 
automatic border B, the detection error was calculated as the mean probability of misclassification over the pixels inside the 
border, i.e. ^fp^pp X 100%. 

Iyatomi et al. [21] modified Guillod et a/.'s approach by combining multiple manual borders that correspond to each lesion 
into one using the majority vote rule. The automatic borders were then compared against these combined ground-truth images. 
Cclcbi et al. [121 compared each automatic border against multiple manual borders independently 

Unfortunately, the above-mentioned methods do not accurately capture the variations in the manual borders. For example, 
according to Guillod et a/.'s measure an automated border that is entirely inside the manual borders would get a very low 
error. Iyatomi et a/.'s method discounts the variation in the manual borders by simple majority voting, while Celebi et a/.'s 
approach does not produce a scalar error value, which makes comparisons more difficult. 

Recently, Celebi et al. [53] proposed the use of an objective measure, the Normalized Probabilistic Rand Index [35J, which 
takes into account the variations in the manual borders. They demonstrated that the differences between four of the evaluated 
border detection methods were in fact smaller than those predicted by the commonly used XOR measure. Since the formulation 
of this measure is involved, the interested reader is referred to and [25] . 
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None of the above measures quantify the effect of border detection error upon the accuracy of the classifier. Loss of 
classification accuracy due to automatic border error can be measured as the difference of the classification accuracy using the 
manual borders and that using the automatic borders. 

6 Comparisons and Discussion 

Table [5] compares some of the recent border detection methods based on their operation mode (automatic vs. semi-automatic), 
color space, preprocessing steps, and segmentation method. Note that only those methods that are adequately described in the 
literature are included and the postprocessing steps are omitted since they are often not reported. The following observations 
are in order: 

• 16/18 methods are automated, reflecting the need for an unsupervised decision support system for clinical use. 

• 11/18 methods operate on multiple color channels. 

• 12/18 methods use a smoothing filter (the hair removal method described in |15j also employs image smoothing). 

• Clustering is the most popular segmentation method, which is probably due to the availability of robust algorithms. 

Table [3] compares the border detection methods based on their evaluation methodology: the number of human experts who 
determined the manual borders, the number of images used in the evaluations (and the diagnostic distribution of these images 
if available), the number of automated methods used in the comparisons, and the measure used to quantify the border detection 
error. It can be seen that: 

• 9/18 studies rely on borders determined by a single dermatologist. 

• Only 5/18 studies report the diagnostic distribution of their test images. This information is valuable given that not 
every diagnostic class is equally challenging from a border detection perspective. For example, it is often more difficult 
to detect the borders of melanomas and dysplastic nevi due to their irregular and fuzzy (hazy) border structure. 

• 8/18 studies do not compare their results to those of any other automated method. This is partly due to the unavailability 
of public border detection software, as well as the non-existence of a public dermoscopy image database. 

• Recent studies used objective measures to determine the validity of their results, whereas earlier studies relied on visual 
assessment. XOR measure is the most popular objective error function despite the fact that it is not trivial to extend 
this measure to capture the variations in multiple manual borders. 

We believe that in a systematic border detection study: 

1. The image acquisition procedure should be described in sufficient detail. 

2. The test image set should be selected randomly from a large and diverse image database. 

3. The test image set should be large enough to ensure statistically valid conclusions. 

4. The diagnostic distribution of the test image set should be stated. 
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Table 2: Characteristics of border detection methods (nr: not reported, KL{C}: KL transform of the C color space) 



Ref. 


Year 


Mode 


Color Space (# Channels) 


Preprocessing Steps 


Segmentation Method 


M 


2008 


Auto 


RGB (3) 


Marginal MF, Black frame removal 


Region- Based 


M 


2008 


Auto 


L*a*b* (3) 


nr 


Clustering 


E 


2008 


Auto 


IHP{RGB} (1) 


Contrast enhancement, Hair removal |15j 


Clustering 


ES 


2007 


Semi 


Luminance (1) 


GF 


Active Contours 


M 


2007 


Auto 


RGB (3) 


nr 


Active Contours 


M 


2007 


Auto 


RGB (3) 


Marginal MF, Approx. Lesion Localization 


Region-Based 


EU 


2006 


Auto 


B / RGB (1) 


GF 


Thresholding 


ED] 


2006 


Auto 


RGB (3) 


nr 


Clustering 


E3 


2005 


Auto 


Luminance (1) 


GF, Auto. Snake Initialization 


Active Contours 


EH] 


2003 


Auto 


L*u*v* (3) 


Marginal MF 


Clustering 


E3 


2002 


Auto 


Kh{L*a*b*} (2) 


GF 


Clustering 


H3 


2001 


Auto 


KL{RGB} (1) 


Marginal MF 


Thresholding 


EH 


2000 


Semi 


L*a*b* (3) 


nr 


Active Contours 


m 


2000 


Auto 


I / HSI (1) 


MF 


Soft Computing 


m 


1999 


Auto 


KL{L*u*v*} (2) 


MF 


Clustering 


m 


1999 


Auto 


L*a*b* (3) 


Vector ADF, Hair removal 


Morphological 


ED 


1998 


Auto 


RGB (3) 


nr 


Region-Based 


EH 


1998 


Auto 


KL{RGB}(1) 


nr 


Model-Based 



5. Algorithms with reasonable computational requirements should be used. 

6. The results should be evaluated using borders determined by multiple dermatologists. 

7. The results should be compared to those of published border detection methods. 

8. The border detection procedure should be described in sufficient detail. 

Note that all of the abovemcntioned criteria except for (5) and (8) can be satisfied by using a public dermoscopy image set. 
Therefore, the creation of such a benchmark database should be prioritized in order to improve the quality of future border 
detection studies. 

We must recognize the limitations of automated systems at present in comparison to experienced dermoscopy practitioners, 
who have the unparalleled ability to correctly identify borders in several situations. First, 'collision tumors', contiguity of 
lesions of more than one type, are fairly common. The most common of these collisions and the situation most likely to cause 
error in automatic diagnosis is the collision between a malignancy and a benign lesion, often a lentigo, which is common in 
sun-damaged skin. Dermatologists arc better able to define what is normal for a given patient's background skin, even when 
the background skin includes lentigines that can falsely enlarge the border. Second, some significant melanoma features may 
be lost without higher level knowledge that enables the inclusion of such features within the lesion. Scar-like regression, an area 
of pallor that is a significant feature for melanoma in-situ, is excluded from the lesion if it is close to the border in all border 
detection methods that we have examined. Third, dermatologists vary the borders according to the diagnosis. Dermatologist 
borders include the halo in halo nevi, and the pale rim of basal cell carcinomas, but exclude the surrounding reactive erythema 
of irritated lesions. It is likely that some higher-level knowledge will need to be included in automated border detection methods 
to accomplish the ultimate purpose of these computer systems: to achieve higher diagnostic accuracy. 
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Table 3: Evaluation of border detection methods (b: benign, m: melanoma) 



Ref. 


# Experts 


# Images (Distr.) 


# Comp. 


Error Measure (Value) 


ma 


3 


90 (65 b / 25 m) 


4 


XOR (10.63%) 




1 


67 





XOR (14.63%) 


m 


1 


100 (70 b / 30 m) 


3 


XOR (2.73%) 


m 


1 


50 


2 


Error prob. (16%) 


M 


1 


50 


2 


Error prob. (21%) 


M 


2 


100 (70 b / 30 m) 


3 


XOR (12.02%) 


M 


5 


319 (244 b / 75 m) 


1 


Prcc. (94.1%) & Rec. (95.2%) 


20 


nr 


117 


3 


Sens. (95%) & Spec. (96%) 




2 


100 (70 b / 30 m) 


1 


XOR (15.59%) 


M 





nr 





nr 


29 





600 





Visual 


30 





nr 





nr 


M 


5 


30 





Visual 


m 


1 


30 





Visual 


CD 


1 


400 





Visual 


m 


1 


300 





Visual 


M 


1 


57 


5 


XOR (36.50%) 


m 


1 


57 


5 


XOR (24.71%) 



Acknowledgments 

This publication was made possible by grants from The Louisiana Board of Regents (LEQSF2008-11-RD-A-12) and The 
National Institutes of Health (SBIR #2R44 CA-101639-02A2). Its contents are solely the responsibility of the authors and do 
not necessarily represent the official views of The Louisiana Board of Regents or the NIH. 



References 



[1] Jcmal A, Siegel R, Ward E et al. Cancer Statistics, 2008. CA: A Cancer Journal for Clinicians 2008; 58(2): 71-96, 2008. 

[2] Argenziano G, Soyer HP, De Giorgi V et al. Dermoscopy: A Tutorial. Milan, Italy: EDRA Medical Publishing & New 
Media, 2002. 

[3] Menzies SW, Crotty KA, Ingwar C, McCarthy WH. An Atlas of Surface Microscopy of Pigmented Skin Lesions: Der- 
moscopy. Sydney, Australia: McGraw-Hill, 2003. 

[4] Steiner K, Binder M, Schemper M et al. Statistical Evaluation of Epiluminescence Dermoscopy Criteria for Melanocytic 
Pigmented Lesions. Journal of American Academy of Dermatology 1993; 29(4): 581-588. 

[5] Binder M, Schwarz M, Winkler A et al. Epiluminescence Microscopy. A Useful Tool for the Diagnosis of Pigmented Skin 
Lesions for Formally Trained Dermatologists. Archives of Dermatology 1995; 131(3): 286-291. 

[6] Fleming MG, Stcger C, Zhang J et al. Techniques for a Structural Analysis of Dcrmatoscopic Imagery. Computerized 
Medical Imaging and Graphics 1998; 22(5): 375-389. 

[7] Celebi ME, Kingravi HA, Uddin B et al. A Methodological Approach to the Classification of Dermoscopy Images. Com- 
puterized Medical Imaging and Graphics 2007; 31(6): 362-373. 

[8] Pratt WK. Digital Image Processing: PIKS Inside. Hoboken, NJ: John Wiley & Sons, 2007. 

[9] Delgado D, Butakoff C, Ersboll BK, Stoecker WV. Independent Histogram Pursuit for Segmentation of Skin Lesions. 
IEEE Trans, on Biomedical Engineering 2008; 55(1): 157-161. 

[10] Celebi ME, Kingravi HA, Aslandogan YA. Nonlinear Vector Filtering for Impulsive Noise Removal from Color Images. 
Journal of Electronic Imaging 2007; 16(3): 033008 (21 pages). 



Computerized Medical Imaging and Graphics, 33(2): 148-153, 2009 



[11] Schmid P. Segmentation of Digitized Dcrmatoscopic Images by Two-Dimensional Color Clustering. IEEE Trans, on Medical 
Imaging 1999; 18(2): 164-171. 

[12] Celebi ME, Kingravi HA, Iyatomi H et al. Border Detection in Dermoscopy Images Using Statistical Region Merging. Skin 
Research and Technology 2008; 14(3): 347-353. 

[13] Geusebroek J-M, Smoulders AWM, van de Weijer J. Fast Anisotropic Gauss Filtering. IEEE Trans, on Image Processing 
2003; 12(8): 938-943. 

[14] Perreault S, Hebert P. Median Filtering in Constant Time. IEEE Trans, on Image Processing 2007, 16(9): 2389-2394. 

[15] Lee TK, Ng V, Gallagher R. et al. Dullrazor: A Software Approach to Hair Removal from Images. Computer in Biology 
and Medicine 1997; 27(6): 533-543. 

[16] Zhou H, Chen M, Gass R et al. Feature-Preserving Artifact Removal from Dermoscopy Images. Proc. of the SPIE Medical 
Imaging 2008 Conf. 6914: 69141B-69141B-9. 

[17] Wighton P, Lee TK, Atkins MS Dcrmascopic Hair Disocclusion Using Inpainting. Proc. of the SPIE Medical Imaging 2008 
Conf. 6914: 691427-691427-8. 

[18] Sonka M, Hlavac V, Boyle R. Image Processing, Analysis, and Machine Vision. Cengage-Engineering, 2007. 

[19] Celebi ME, Aslandogan YA, Stoecker WV et al. Unsupervised Border Detection in Dermoscopy Images. Skin Research 
and Technology 2007; 13(4): 454-462. 

[20] Melli R, Grana C, Cucchiara R. Comparison of Color Clustering Algorithms for Segmentation of Dcrmatological Images. 
Proc. of the SPIE Medical Imaging 2006 Conf. pp. 3S1-9. 

[21] Iyatomi H, Oka H, Saito M et al. Quantitative Assessment of Tumor Extraction from Dermoscopy Images and Evaluation 
of Computer-based Extraction Methods for Automatic Melanoma Diagnostic System. Melanoma Research 2006; 16(2): 
183-190. 

[22] Hancc GA, Umbaugh SE, Moss RH, Stoecker WV. Unsupervised Color Image Segmentation with Application to Skin 
Tumor Borders. IEEE Engineering in Medicine and Biology 1996; 15(1): 104-111. 

[23] Guillod J, Schmid-Saugeon P, Guggisberg D et al. Validation of Segmentation Techniques for Digital Dermoscopy. Skin 
Research and Technology 2002; 8(4): 240-249. 

[24] Celebi ME, Iyatomi H, Schaefer G et al. An Improved Objective Evaluation Measure for Border Detection in Dermoscopy 
Images. Skin Research and Technology 2009; 15(4): 444-450. 

[25] Unnikrishnan R, Pantofaru C, Hcbcrt M. Toward Objective Evaluation of Image Segmentation Algorithms. IEEE Trans, 
on Pattern Analysis and Machine Intelligence 2007; 29(6): 929-944. 

[26] Mendonca T, Marcal ARS, Vieira A et al. Comparison of Segmentation Methods for Automatic Diagnosis of Dermoscopy 
Images. Proc. of the 29th IEEE EMBS Annual Int. Conf. 1: 6572-6575. 

[27] Erkol B, Moss RH, Stanley RJ, Stoecker WV, Hvatum E. Automatic Lesion Boundary Detection in Dermoscopy Images 
Using Gradient Vector Flow Snakes. Skin Research and Technology 2005; 11(1): 17-26. 

[28] Galda H, Murao H, Tamaki H, Kitamura S. Skin Image Segmentation Using a Self-Organizing Map and Genetic Algorithms. 
Trans, of the Institute of Electrical Engineers of Japan - Part C 2003; 123(11): 2056-2062. 

[29] Cucchiara R, Grana C, Scidenari S, Pcllacani G. Exploiting Color and Topological Features for Region Segmentation with 
Recursive Fuzzy c-means. Machine Graphics and Vision 2002; 11(2/3): 169-182. 

[30] Hintz-Madsen M, Hansen LK, Larsen J, Drzewiecki K. A Probabilistic Neural Network Framework for the Detection of 
Malignant Melanoma. Artificial Neural Networks in Cancer Diagnosis, Prognosis and Patient Management, pp. 141-183, 
2001. 

[31] Haeghen YV, Naeyaert JM, Lemahieu I. Development of a Dermatological Workstation: Preliminary Results on Lesion 
Segmentation in CIE L*A*B* Color Space. Proc. of the Int. Conf. on Color in Graphics and Image Processing. 

[32] Donadey T, Serruys C, Giron A et al. Boundary Detection of Black Skin Tumors Using an Adaptive Radial-based Approach. 
Proc. of the SPIE Medical Imaging 2000 Conf. 3379: 810-816. 

[33] Schmid P. Lesion Detection in Dcrmatoscopic Images Using Anisotropic Diffusion and Morphological Flooding. Proc. of 
the IEEE ICIP 1999 Conf. 3: 449-453. 

[34] Gao J, Zhang J, Fleming MG et al. Segmentation of Dermatoscopic Images by Stabilized Inverse Diffusion Equations. 
Proc. of the IEEE ICIP 1998 Conf., 3: 823-827. 



